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EFFECT OF ICE CONTAMINATION ON LIQUID -NITROGEN 
DROPS IN FILM BOILING 
by G. J. Schoessow,* C. E. Chmielewski, ** 
and K, J. Baumeister 

NASA Lewis Research Center 
Cleveland, Ohio 

ABSTRACT 

Previously reported vaporization time data of liquid nitrogen drops 
in film boiling on a flat plate are about 30 percent shorter than predicted 
from standard laminar film boiling theory. This theory, however, had 
been found to successfully correlate the data for conventional fluids such 
as water, ethanol, benzene, or carbon tetrachloride. This paper presents 
experimental evidence that some of the discrepancy for cryogenic fluids 
results from ice contamination due to condensation. The data indicate 
a fairly linear decrease in droplet evaporation time with the diameter 
of the ice crystal residue. After correcting the raw data for ice contam- 
ination along with convection, a comparison of theory with experiment 
shows good agreement, 

INTRODUCTION 

Leidenfrost film boiling data of discrete masses of a cryogenic fluid 
were first reported by Keshock and Bell (ref, 1 and 2). They attempted 
to correlate their cryogenic data using the equations developed in refer- 
ence 3 from standard laminar film boiling theory. This theory had been 
found to successfully correlate data for conventional fluids such as water, 
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ethanol, benzene, or carbon tetrachloride. However, Keshock and Bell 
found the measured drop vaporization times were as much as 30 percent 
shorter than predicted by the equations developed in reference 3, even 
after the contributions of heat transfer by radiation and free convection 
to the droplet's upper surface were considered. To account for this 
discrepancy, the effect of ice contamination and mass transfer from 
the upper surface of the drop were experimentally investigated in the 
present paper. In particular, ice formation in the drop from contam- 
ination was found to significantly affect the vaporization times of liquid 
nitrogen drops. 

INITIAL EXPERIMENTS 

Keshock and Bell (ref. 2) suggested that the difference between 
theory and experiment, particularity at large drop volumes, results 
from vapor break through beneath the drop, similar to that observed 
with pool film boiling. Later, Baumeister, Keshock and Pucci (ref. 4) 
used this hypothesis to develop a correlation based on a heat transfer 
coefficient which shifted from a droplet coefficient to a pool boiling 
coefficient when the drop volume was sufficiently large so that vapor 
break through occurred. This technique produces some improvement 
between theory and experiment for large nitrogen drops (see ref. 4, 
fig. 1); however, the measured evaporation times were still lower than 
those predicted by theory. 

As part of this paper, a 5 cm diameter nitrogen cooled shield was 
placed above the tops of the drops (within 0. 5 cm) to reduce (eliminate) 
radiation, convection and mass transfer from the upper drop surface, 
much the same as a nitrogen shield used in dewars for storage and 
transport of liquid hydrogen or helium. For liquid nitrogen drops on 
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a plate at 20° C, no significant change in the total vaporization times were 
seen, ‘in comparison to the vaporization times without the shield. However, 
in special cases where radiation and convection could be important (moving 
drop) Keshock and Bell report (ref. 1) that even after correction for radia- 
tion and free convection are made, the experimental results are still sub- 
stantially lower than theoretical predictions. 

During these initial experiments, a problem was encountered with 
the absorption of water into the drop and the formation of a slight ice 
residue. Keshock, in his very exacting work (ref. 4), dehumidified the 
droplet environment with dessicants and still had problems with ice con- 
tamination in the droplets. There was therefore, the possibility that 
previous data taken for film boiling studies of nitrogen were effected by 
contamination. The experiments performed herein assess the magnitude 
of this effect on the vaporization times. 

EXPERIMENTAL APPARATUS AND PROCEDURE 

Leidenfrost boiling of drops was performed on a flat stainless steel 
plate (slightly concaved at center to prevent drop motion) and on a vertically 
mounted aluminum wheel (ref. 5) which could be spun to simulate the ef- 
fect of drop velocity on Leidenfrost boiling. Thermocouples were used to 
monitor the surface temperature. The nitrogen was delivered by use of 
two metal cups with insulated handles to reduce heat loss. The cups 
were calibrated by a weighing procedure with distilled water (see ref, 6). 

The diffusive and convection effects by different cover gases appear 
to be quite small (ref, 7). The initial experiments indicated that the 
vaporization times of liquid nitrogen drops in air agreed with Keshock’s 
tlvXlh. (vetv U for a dry nitrogen atmosphere. Because of these observa- 
tions ^ liquid nitrogen drop behavior was investigated in an ambient air 
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atmosphere since this would simplify the procedures and equipment re- 
quirements and allow close observation of the evaporation process, The 
nitrogen was simply delivered to the heat transfer surface using the de- 
livery cup and a timer was started. Extreme care was taken to prevent 
drop fragmentation during delivery. 

During these experiments, the size of the ice residue increased with 
subsequent data runs. However, when a new container of nitrogen was 
obtained, the ice residue size abruptly dropped. Careful observation of 
the liquid nitrogen handling and storage equipment indicated a possible 
source of the contamination was from frost formation on the storage flash 
spout and much of the handling equipment. By air drying all apparatus 
and cleaning the equipment periodically during an experiment, the amount 
of ice crystal impurity present could be held at a nearly constant low 
value but could not be completely eliminated. Some of the contamination 
was apparently coming from storage and possibly the manufacturing process 
itself. 

To investigate the effect of ice contamination on evaporation time, 
frost was deliberately introduced into storage containers of nitrogen 
which were used in an experiment. In the following experiments, all the 
parameters were held constant except the contamination volume. To 
measure the residue ice volume, a micrometer was placed above the flat 
plate approximately 3/8 inch above the droplet. Since the ice mass re- 
maining did not melt for several seconds after the vapor film collapsed, 
there was time to measure (estimate) the diameter with the micrometer. 

The total vaporization time could then be related to the diameter of the 
ice residue, see figure 1, Next, earlier experiments were repeated taking 
care to measure the diameter of the ice residue. 


5 


RESULTS 

The data indicated a fairly linear decrease in evaporation time with 
diameter of the ice residue over the range where the residue appears 
to be fairly spherical, as shown in figure 1. Next, the measured evap- 
oration times of the earlier experiments were corrected for ice contamina- 
tion using a linear least squares data fit. The complete procedure for 
correcting the raw data to account for ice formation is given in Appendix C 
of ref. 6. A comparison of the theory to the corrected vaporization times 
are shown in figure 2. The data were also corrected for free convection 
using the Keshock and Bell procedure (refs, 1 and 2). The data fall into 
the extended drop range shown in table I. 

The theoretical relationship for the vaporization time, shown in fig- 
ure 2, was developed in ref. 8 by first solving the Laplace capillary 
equation for the drop thickness and radius as a function of volume. A closed 
form analytical solution did not exist over the complete range of drop volume; 
however, by asymptotic assproximations, closed form solutions were de- 
veloped for three separate volume ranges - small, intermediate, large. 

Next, these solutions for drop shape were combined with the simultaneous 
solution of the continuity, momentum, and energy equations in the laminar 
vapor layer beneath the drop to obtain the heat transfer coefficient between 
the supporting plate and the drop. Finally, a time dependent energy balance 
on the droplet lead to expression for the total vaporization time of the 
drops. Because the drop shape was broken into three separate ranges, 
three separate solutions for the vaporization time are required. These 
final expressions are shown in table I, 

As seen in figure 2, good agreement exists between experiment and 
theory. The correction for ice formation and convection is approximately 
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20 .percent over a range of drop velocities of 0 to 3. 7 m/sec. For larger 
drop volumes, the correction can be expected to be even larger, as pointed 
but in reference 2. 

CONCLUSIONS 

When convection, bubble break through, and ice contamination correc- 
tions are applied to vaporization time data of cryogenic nitrogen drops in 
Leidenfrost boiling, predictions of the drop vaporization times based on 
standard laminar film boiling theory are in reasonable agreement with 
experiment. 

NOTATION 

2 constants, see table I 

Cp specific heat at constant pressure of vapor 
F dimensionless function, defined in table I 

g coefficient of gravity 

g gravitational constant (conversion constant between mass and 

force units) 

k thermal conductivity of vapor beneath drop 

I average drop thickness, defined table I 

Z* dimensionless average drop thickness, see table I 

P Prandtl number of vapor, C ju/k 

X y 

AT temperature difference, T - T„ 

w s 

T„ wall temperature 

T„ saturation temperature 

O 

t time 

t* dimensionless time, defined in table I 

U velocity of plate relative to drop 

V volume of drop 
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V* dimensionless volume of drop, see table I 

r see table I for definition 

X latent heat of vaporization / 

/ 7 C AT 

X* modified heat of vaporization A* = X i + — -° — 

\ 20 A 

Ij. viscosity of vapor 

pQ. vapor density 

o 

pj^ liquid density 

latent to sensible heat ratio, X*/C AT 

It’ 
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TABLE I. - SUMMARY TABLE FOR MOVING DROP RESULTS 


spheroid 
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’777777 


<tC2S> 



Drop thickness, 

Vaporization time. 






l^kW'^^(PL-PB)Pg‘iT2/ 

I* = 0. 83 

t* = 1.21 V+®/^^ 

I* = 0, 8 

t* = 2.23 - 0.37 

1 * ^ 1.85 

t» = 4. 52 - 5 


Pseudo velocity dependent volume, V 



V+ = V* (r(V) + 


4,52 V»‘ 
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Figure 1. - Effect of ice crystal contamination on evaporation lime of 0.47 ml stationary nitrogen droplets on a 
flat plate. 
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Figure 2. - Dimensionless total vaporization correlallon of nitrogen data with 
Ice and free convection corrections. 


